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Abstract 

A preliminary determination of the first moment of the hadronic mass distribution (M^ — Tn 2 £>) in 
semileptonic B decays has been obtained as a function of the minimum lepton momentum, ranging 
from 0.9 to 1.6 GeV/c. The measurement is based on a new technique involving BB events in which 
one fully reconstructed B meson decays hadronically and the recoiling B decays semileptonically. 
The mass of the hadrons in the semileptonic decay is determined from a kinematic fit to the 
whole event. For different minimum lepton momenta, the mass distribution is decomposed into 
contributions from various charm resonant states and a non-resonant contribution, allowing for 
the determination of the first moment. From these moments the Heavy Quark Effective Theory 
(HQET) parameters Ai and A can be derived. For lepton momenta in the B rest frame above 
1.5GeV/c, we find a first moment that is compatible with existing measurements. However, if we 
extend the measurement to lower values of lepton momenta, the data can only be described by 
Operator Product Expansion calculations if we use significantly different values for A and Ai than 
obtained from earlier measurements based on lepton momentum spectra and the photon spectrum 
in b — > sj transitions. 
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1 Introduction 



The heavy quark limit in QCD has become a very useful tool for relating inclusive B decay proper- 
ties, like the semileptonic branching fraction, to the charged current couplings, |V^,| and \V u i\. These 
inclusive observables can be calculated using expansions in powers of the strong coupling constant 
a s (mb) and in inverse powers of the B meson mass, tub, that include non-perturbative parameters 
A, Ai and A2. From the B* — B mass splitting A2 has been determined to be 0.128 ± 0.010 GeV 2 
and it is expected that A and Ai will be calculated with techniques such as lattice QCD. However, 
they can also be extracted experimentally from inclusive measurements. 

In this paper we report a measurement of the invariant mass Mx of the hadronic system 
recoiling against the charged lepton and the neutrino in semileptonic B decays. In the framework 
of HQET, calculations of the moment (M x — fn 2 D ) have been carried out |Q] using an Operator 
Product Expansion (OPE) in inverse powers of the B meson mass up to order a 2 /?o and l/ro|. Here 
mo = (mD+3mo«)/4 = 1.975 GeV/c 2 is the spin-averaged D meson mass and (3q = (33 — 2nj)/3 = 
25/3 is the leading order QCD (3 function. In these expansions the non-perturbative parameters A 
(0(1/ tub)), Ai and A2 {0(l/m 2 B )) appear. A measurement over a large region of phase space can 
test the validity of the assumptions made in the OPE, such as parton-hadron duality. 

Our determination of the hadronic moment in semileptonic B decays exploits new capabilities 
made possible by very large samples of reconstructed B mesons, in contrast to earlier measurements 
by CLEO (||. The data sample for the analysis relies on BB events in which one B meson is fully 
reconstructed and thus the momentum of the semileptonically decaying B is known. This allows us 
to assign the remaining tracks and clusters in the event to the hadronic mass Mx of the semileptonic 
decay. Furthermore, the determination of Mx can be substantially improved by application of a 
two-constraint kinematic fit to the full event. 

The moment {M x —rfijj) of the Mx distribution is measured for lepton momenta (calculated 
in the B rest-frame) as low as 0.9 GeV/c. The wide momentum range improves the sensitivity of 
the measurement to a variety of hadronic final states, including higher resonances of the D and 
non-resonant D^ir combinations. 



2 The BABAR Detector and Data Set 

Our analysis is based on data recorded with the BABAR detector |J at the PEP-II energy-asymmetric 
e + e~ storage ring at SLAC. The detector consists of a five-layer silicon vertex tracker (SVT), a 40- 
layer drift chamber (DCH), a detector of internally reflected Cherenkov light (DIRC), and an elec- 
tromagnetic calorimeter (EMC) assembled from 6580 CsI(Tl) crystals, all embedded in a solenoidal 
magnetic field of 1.5 T and surrounded by an instrumented flux return (IFR). The data sample 
corresponds to an integrated luminosity of 51 fb _1 that was collected at the T(4S) resonance. 

3 Analysis Method 

The basic analysis technique is based on BB events in which one of the B mesons decays hadron- 
ically and is fully reconstructed (B reco ) and the recoiling B (B reco ii) decays semileptonically. The 
requirement of such a decay topology allows the B momentum to be determined, so that transfor- 
mations to the rest frame of the recoiling B meson are possible. It also results in a clean sample of 
BB events, with the flavor of the reconstructed B meson determined. 
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3.1 Selection of Hadronic B Decays, B — > DY 

The reconstruction of B decays is designed to identify samples of decays of the type B — > DY. Here 
.D refers to a charm meson and Y represents a collection of hadrons of total charge ±1, composed 
of n^, K^, Kg and 7r° mesons. The charm meson serves as a "seed" for the selection of candidates. 
We use four different seeds: D + and D* + for B° and D° and D*° for B~ (charge conjugate modes 
are implied throughout this document). Several different decay modes are used, each of which is 
characterized by a signal-to-background ratio dependent on the multiplicity and on the composition 
of the Y system. In events with more than one reconstructed B decay, the decay mode with the 
highest a priori signal-to-background ratio is selected. 

B reco candidates are identified by two kinematic variables, the beam energy-substituted mass, 

Mes = \j^beam ~ P*b ' an d tne ener gy difference in the T(4S) rest frame, AE = E* B — E£ eam . Here 
El eam is the beam energy, and p* B and E B refer to the B momentum and energy in the center-of- 
mass frame. B reco candidates are required to have AE within three standard deviations of zero 
and M E s > 5.27GeV/c 2 . A sideband region in Mes used for background subtraction is defined by 
5.20 < Mes < 5.27GeV/c 2 . An example for the Mes distribution can be found in Figure |]. 

3.2 Selection of Semileptonic Decays, B — > Xtv 

Semileptonic B decays are identified by the presence of one and only one electron or muon above a 
minimum momentum P^ n measured in the rest frame of the B meson recoiling against the B reco . 
Pmin i s var i e d for the different measurements in the range of 0.9 — 1.6GeV/c. This requirement on 
the lepton momentum reduces backgrounds from secondary charm or decays, and from hadrons 
faking leptons. 

To further improve the background rejection, we require that the lepton charge and the flavor 
of the reconstructed B^ be consistent with a prompt semileptonic decay of the B reco n. For a 
reconstructed B° we do not impose this restriction, so as not to lose events due to B° — B° mixing. 

Electrons are identified using a likelihood-based algorithm, combining the track momentum 
with the energy, position, and shape of the shower measured in the EMC, the Cherenkov angle 
and the number of photons measured in the DIRC, and the specific energy loss in the DCH. The 
efficiency of these selection requirements has been measured with radiative Bhabha events as a 
function of the laboratory momentum pi^ and polar angle 0^, and has been corrected for the 
higher multiplicity of BB events using Monte Carlo simulations. The average electron selection 
efficiency is (90 ±2)%. 

Muons are identified using a cut-based selection for minimum ionizing tracks, relying on informa- 
tion from the finely segmented instrumented flux return of the magnet. The number of interaction 
lengths traversed by the track, the spatial width of the observed signals, and the match between 
the IFR hits and the extrapolated charged track are used in the selection. The muon identification 
efficiency has been measured with ^^(j) events and two-photon production of ^ [i~ pairs. The 
average muon selection efficiency is (70 ± 4)%. 

The misidentification probabilities for pions, kaons, and protons have been extracted from 
selected control samples in data. They vary between 0.05% and 0.1% (1.0% and 3%) for the 
electron (muon) selection. 

To further reduce backgrounds we require the total charge of the event to be \Qtot\ = \Q Breco + 
QBrecoiil < 1; and restrict the total missing mass to | A£^ iss | < 1.0 [ GeV/c 2 ] 2 . We explicitly allow for 
a charge imbalance to reduce the dependence on the exact modeling of charged particle tracking in 
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Figure 1: The Mes distribution for selected events with one reconstructed B decay and a lepton 
with P^ in = 0.9GeV/c. 



the Monte Carlo simulation, especially at low momenta, and the production of tracks from photon 
conversions. 

Figure || shows the Mes distribution for events with a B reco candidate and a charged lepton 
candidate from the recoiling B meson. A sideband subtraction is performed in bins of Mx to 
remove the background contribution under the peak at the B mass. For this purpose a sum of 
background (approximated by a function first introduced by the ARGUS Collaboration J|]) and 
signal (described by a function first used by the Crystal Ball experiment is fitted to the Mes 
distribution of the data. For Mes > 5.27GeV/c 2 we find 5819 signal events and 3585 background 
events. The fit to the Mes distribution is shown in Figure |l|. 

3.3 Reconstruction of the Hadronic Mass M x 

The hadron system in the decay B — » Xtv is formed from K^ 1 , ir^ 1 , and photons that are not 
associated with the B reco candidate or identified as leptons. Specifically, we select charged tracks in 
the fiducial volume 0.41 < 6^ < 2.54 rad with a minimum transverse momentum pt > 120MeV/c. 
For photons we require an energy Ei a b > 100 MeV and 0.41 < 6i a b < 2.409 rad. 

To improve the resolution on the measurement of the hadronic mass Mx , we exploit the kine- 
matic constraints of the BB state by performing a kinematic fit to the full event using the measured 
momenta and energies of all particles. The measured four-momentum P'^ of the X system can be 
written as 
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Figure 2: Monte Carlo comparison of the reconstructed and fitted mass Mj: Mean and r.m.s. of 
the distribution AM X = M\ rue - M^ eco (open circles) and AM X = M\ rue - M™ (dots) as a 
function of the measured missing mass squared M 2 
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where P are four-momenta and the indices ch and 7 refer to the selected charged tracks and 
photons, respectively. Depending on particle identification the charged tracks are assigned either 
the or ir^ mass. 

The missing momentum in the event is reconstructed by relating the sum of the four-momenta 
of the colliding beams, Qcm, to the measured four-momenta of the B reco candidate, the X hadrons, 
and the charged lepton, 

Prmss = QCM - P r m eco ~ P% ~ P™ ■ (2) 

The measured invariant mass squared, M^ iss = P m i ss , is an important discriminant on the quality 
of the reconstruction of the total recoil system. Any secondary particles from the decay of the 
hadronic X system that are undetected or poorly measured will impact the measurement of both 
Mx and M^ niss . Likewise any sizable energy loss of the leptons via bremsstrahlung or internal 
radiation will impact the measurement of these two quantities. The effect of initial state radiation 
is rather small, due to the narrow width of the Y(AS) resonance. We exploit the available kinematic 
information from the full event, namely the B reco and the B recoi i candidate, by performing a 2C 
kinematic fit that imposes four-momentum conservation, the equality of the masses of the two B 
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Figure 3: Sideband- subtracted Mx distribution for Pj^ in = 0.9GeV/c (left) and P^in = l-5GeV/c 
(right). The hatched histograms show the contributions from B — > D*£v, B — ► Dlv and B — > 
Xfttv decays as determined by the fit, as well as the background distribution. The white histogram 
represents the sum of all the aforementioned distributions. The histograms extend over the range 
in Mx that was used in the fit to extract {M x — m 2 D ) . 

mesons, M reco u = M reco , and forces M^ iss = M 2 = 0. The fit takes into account event-by-event 
the measurement errors of all individual particles and the measured missing mass. 

Figure |2| shows a comparison for B — * X c lv events between the reconstructed and the kinemat- 
ically fitted mass Mx as a function of the measured M^ iss , in terms of the mean and r.m.s. of the 
distributions for AMx = Mx (true) — Mx ■ 

The constraints of energy and momentum conservation and the equal-mass hypothesis lead not 
only to a significant improvement in the mass resolution of the X system but also provide an almost 
unbiased estimator of the mean and a resolution that is largely independent of M^ niss . 



4 Determination of the Hadronic Mass Moments 

Figure |3| shows the measured Mx distribution of the selected events after sideband subtraction, 
for two different minimum lepton momenta P^in- Several B — > X c lu decays contribute to this 
distribution. The dominant decays are B — > D*lv and B — > Dlv, but there are also contributions 
from decays to higher mass charm states, D** resonances with a mass distribution X^ so peaked near 
2.4 GeV/c 2 , and non-resonant D^ir final states with a broad mass distribution X^f eso extending 
to the kinematic limit. 
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Table 1: Preliminary results for {M x —m 2 D ) for different values of P^ in . The errors are statistical 
and total systematic respectively. The last four columns show separately the four dominant 
contributions to the systematic uncertainty: Xh model dependence, detector resolution, background 
contribution, and sideband (SB) background subtraction. 
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The observed Mx distribution is fitted to the sum of four distributions that we derive from 
Monte Carlo simulations: three contributions for the decays B — > D*lv, B — > Dlv and B — > Xftlv, 
and one for backgrounds. The ratio of the total branching fractions for the production of X^ so 
and Xjf eso is fixed to 1.25 as assumed in the Monte Carlo event generator. The contribution of the 
different components to the Mx distribution varies as a function of the charged lepton momentum. 
The shape and contribution of the background also vary as a function of -P^m an d are fixed to 
the prediction from the Monte Carlo simulation. The dominant background is due to secondary 
semileptonic charm decays. Contributions from lepton (primarily muon) misidentification are much 
smaller. The backgrounds decrease significantly for higher P^ in . 

Detailed Monte Carlo simulations are used to reproduce the impact of the efficiencies on signal 
and background and to derive the shape of the signal and background contributions to the Mx 
distribution. These simulations take into account variations of the detector performance and the 
beam background. The decays B — > Dlv and B — > D**lv are simulated using the ISGW2 model g, 
whereas for the decays B — > D*lv a parameterization of HQET-derived form factors is used. 
The non-resonant final states such as Dit and D*n are modeled according to the prescription of 
Goity and Roberts ||. 

A binned x 2 ~fit to the Mx distribution is performed to determine the relative normalization 
of the three signal contributions, fo*, fo-, and fx H - Figure |3] shows the fit results for a minimum 
lepton momentum of -P^ in = 0.9GeV/c and -P^ in = 1.5GeV/c. At higher lepton momenta, the 
contribution from the higher mass states is reduced, specifically the non-resonant charm states 
contribute very little above 1.5GeV/c. For P^ in = 0.9GeV/c decays to D*, D, and higher mass 
final states X r £ so and X% reso contribute approximately 50%, 20%, 20% and 10% of the total M x 
distribution. These fractions change to 60%, 20%, 15% and 5% for P^ in = 1.5GeV/c. Thus the 
average (Mx) value decreases for higher lepton momenta. 
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Figure 4: Measured values for {M'j < _—rn 2 D ) for different minimum lepton momenta, P^in- The inner 
error bars indicate the statistical uncertainty, the total error bars show the quadratic sum of the 
statistical and systematic errors. The errors of the individual measurements are highly correlated. 
The CLEO measurement jj^J is shown for comparison. The lower line shows the predicted variation 
of the moments for the parameters Ai = —0.17 GeV 2 and A = 0.35 GeV ^j, that were obtained from 
the single measurement at P^in = 1-5 GeV/c. For comparison, the dashed curve shows an example 
of an OPE parameterization which follows the P^m dependence of the measured moments. 



14 



QQ 



6 
5 
4 
3 
2 
1 





I 



T 



T 

1 



I 



I 

1 



f > f f f f n 



PDG 



D* 



D 



n x 



J I I I L 



0.8 0.9 1 1.1 1.2 1.3 1.4 1.5 1.6 1.7 

P' min (GeV/c) 



^ 2.5 
m 

2 
1.5 
1 

0.5 




X H (resonant) 

I T I 



X H (non — resonant) 



I 



l 



J I I I L 



0.8 0.9 1 1.1 1.2 1.3 1.4 1.5 1.6 1.7 

P* min (GeV/c) 



Figure 5: The upper figure shows the measured branching fractions of D* , D and Xh as a func- 
tion of Pmin that are obtained from the Mx fits after acceptance correction (points), and the world 
average J7ZJ/ branching fractions for D* and D (yellow bands). The lower figure shows the extracted 
branching fractions obtained from a fit where the resonant and non-resonant components are al- 
lowed to float independently. All branching fractions are obtained by assuming a total semileptonic 
branching fraction of 1 0. 87%. 

The first moment of M x (or any other moment) can be extracted from the fitted relative 
contributions of the different decay modes, /d*, /d, and fx H - These fractions are normalized to 
the B — ► D*lu branching fraction, and as a result many uncertainties in the detection efficiencies 
and branching fractions cancel. Taking into account the true particle masses (the D and D* masses 
are basically 5 functions, and the mean of the Xh contribution is taken from generated events) the 
first moments are calculated according to the following expression 



m 2 D ) 



Id* ■ {Ml, - m 2 D ) + f D • {Ml - m z D ) + f Xa ■ (Ml - ml). 



™2 



™2 



(3) 



The results obtained for (M x — m 2 D ) from the x 2 fit are summarized in Table |T| and Figure |3| as 
a function of the minimum lepton momentum P^ n , ranging from 0.9 — 1.6GeV/c. The errors for 
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different P^ in are highly correlated because all events selected for a certain P^ in are contained in 
samples with lower P min - The CLEO Collaboration has also measured the first hadronic moment 
for P^ in = 1.5 GeV/c. Their result of (M x -m 2 D ) = 0.251 ±0.066 [GeV/c 2 ] 2 § is compatible with 
the measurement presented here. 



5 Studies of Systematic Uncertainties 

Extensive studies have been performed to assess the systematic uncertainties in the measurement 
of the first hadronic moment as a function of P m i n - These studies involve both changes in the event 
selection and variations of the various corrections for efficiencies and resolution, and comparisons 
of data with Monte Carlo simulations for various control samples. The systematic errors are listed 
in Table |. 

With the current data set the total error on the moments is statistically dominated. One of 
the most important systematic errors comes from the uncertainty in the model for the higher mass 
states Xh- The resonant (ISGW2) and non-resonant (Goity-Roberts) decays are each modeled for 
four decay modes. The systematic error is evaluated by alternately turning on and off all individual 
components of the two models and choosing all possible combinations of decay modes for each of 
the two models. In this way not only the shape of the resonant and non-resonant component is 
varied but also the relative branching fractions between the different high-mass states and their true 
underlying masses are altered. We assign the r.m.s. deviation of these variations as the systematic 
error due to the Xh modeling. This error has the largest contribution for low P^ in and decreases 
for higher lepton momenta because the high-mass final states are significantly reduced. 

Another systematic error arises from differences in the track and photon resolution between 
Monte Carlo simulations and control samples in data. The effect on the moments is estimated by 
varying the resolution and changing the bin width of the Mx distribution. 

The third important systematic error is due to the uncertainties in the background contribution 
from secondary semileptonic decays of charm particles. This background was varied by up to ±30% 
in steps of 5% and the x 2 -P r °bability weighted r.m.s. deviation of these variations is taken as the 
systematic error. This error also decreases for higher P m i n - The Mx distribution obtained from 
wrong-sign leptons in data and Monte Carlo simulation is found to be in good agreement. The 
systematic uncertainty due to the sideband subtraction is evaluated by varying the lower border of 
the signal window from 5.265 — 5.275 GeV/c 2 and reducing the sideband region from the default size 
of 5.20 — 5.27 GeV/c 2 to 5.24 — 5.27 GeV/c 2 . This addresses possible Mes dependent changes in 
the shape of the Mx distribution obtained from the sideband as well as uncertainties regarding the 
ratio of background in the signal region to background in the sideband region. The full variation 
of the moments was taken as the systematic error. 

The impact of changing the Mx range for the fit was also studied and the variation of the 
measured moments is taken as a systematic error. Its size is almost P^ in independent and is 
included in the error quoted for detector systematics. 

Several cross checks were carried out. In particular, we performed a model independent, direct 
measurement of Mx for the resonant and non-resonant states by subtracting the known D, D* 
contributions and background from the Mx distribution of the data. The branching fractions for 
D and D* were taken from the PDG [10] and the P min dependent acceptance of both decay modes 



as well as the background were taken from the Monte Carlo simulation. Figure |6j shows the Mx 
distribution before and after subtraction for the lowest P m i n of 0.9 GeV/c. 
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Figure 6: The upper figure shows the sideband- subtracted Mx distribution along with the ex- 
pected D, D* and background components for P^ in = 0.9 GeV/c. The lower figure shows the Mx 
distribution of the high-mass final states obtained after subtracting the D, D* , and background com- 
ponents. The filled histograms correspond to the Monte Carlo expectations for the resonant (grey) 
and non-resonant (black) components. ^ 
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Figure 7: A comparison of the results for the first hadronic moment obtained from pairwise 
independent subsamples. 

A significant contribution to the Mx distribution from high-mass final states can be observed 
and is found to be in good agreement with the Monte Carlo simulation. This measurement of the 
high-mass states is not corrected for detector resolution effects. However, through the use of the 
two-constraint kinematic fit, the mean is believed to be an almost unbiased estimator of its true 
underlying value (see Figure ^) . The mean of the measured mass distribution is in agreement with 
the mean {Mx H ) of the Monte Carlo distribution used in the standard extraction of the moments 
(see Eq. |). 

We also used the relative contributions fi to determine branching fractions by correcting for 
acceptance for each P^ in value and assuming a total semileptonic branching fraction of 10.87% [11|. 
The dependence of these branching fractions on -P^j n is shown in Figure |5[ We obtain branching 
fractions for D and D* that are both compatible with previous measurements [[h| and stable as a 
function of P^im although of course correlated. 

In order to further explore the non-resonant contributions to the Mx distribution we split the 
combined shape used for the Xh system into its resonant and non-resonant components and repeat 
the fits at different -P^m with four independent contributions. The branching fractions obtained 
by this means are also shown in Figure |5[ 

With the current data sample the sensitivity to the four independent contributions is not suf- 
ficient at large -P^m where the non-resonant contribution is suppressed. We find stable results for 
the resonant and non-resonant branching fractions and the ratio between the two is compatible 
with that used in the Monte Carlo simulation. 
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Other tests include a comparison of the moment measurements for different subsamples of the 
data: versus B°, decays with electrons versus muons, and events with total charge Qtot = 
and | Qtot | = 1. A comparison of these results is shown in Figure [7|. 

The effect of variation of the signal-to-background ratio in the B reco samples was also studied. 
Likewise, the requirement on |M^ is J was varied. All results were found to be in good agreement 
with the standard extraction of the moments for the full data sample. 



6 Conclusions 

We have performed a preliminary measurement of the first moment (Mj^ — m 2 D ) of the hadronic 
mass distribution in semileptonic B decays and extend the measurements to lepton momenta as 
low as 0.9 GeV/c. As a result, these measurements are sensitive to higher-mass charm states, in the 
form of both resonant and non-resonant contributions. The full reconstruction of one of the two B 
mesons and the kinematic fit to the full event have resulted in measurements that are limited by 
statistical uncertainties. For lepton momenta above 1.5 GeV/c, the measured moment agrees well 
with an earlier result reported by the CLEO Collaboration [g]. 

Using the OPE, the measurements of (M^ — rn 2 £>) can be translated into values for Ai and A. 
The coefficients used in the expansions are obtained from calculations by Falk and Luke []l], [12| . 
The observed dependence of the moments on the minimum lepton momentum can be reproduced 
by the OPE, as long as we adjust the non-perturbative parameters (see Figure |] for an example). 

If, on the other hand, we insert into the OPE calculations recent independent measurements 
of A, the data cannot be reproduced. The CLEO Collaboration reported a measurement of A = 
0.35 ± 0.08 ± 0.10 GeV ||, derived from the first moment of the photon energy spectrum above 
2.0 GeV in b — > sj transitions. CLEO also recently reported a measurement of A from the lepton 
energy spectrum in semileptonic B decays, for P^ in = 1.5 GeV/c, which yielded a comparable value 



For P^ in = 1.5 GeV/c and A = 0.35 ± 0.13 GeV, we find 

\ 1 = -0.17 ± 0.06 ± 0.07 GeV 2 , (4) 

where the first error is experimental and the second theoretical. As stated above, this result is also 
in good agreement with the CLEO measurement of Ai = —0.236 ± 0.071 ± 0.078 GeV 2 . However, 
when we use the same value of A and the moments measured over a wider lepton momentum range, 
we find values for Ai that are not compatible with each other. Specifically, we obtain 

*lUv =o. 9G ev/ c - A iU. =, 5G ev/ c = - 05 GeV 2 -(-0.17 GeV 2 ) = 0.22 ± 0.04 ± 0.05 GeV 2 . (5) 

The first error is statistical and takes into account the correlation between the individual moment 
measurements. The second error is systematic and is dominated by the uncertainty on A. 

Another way to look at this is to use A x = -0.17 GeV 2 and A = 0.35 GeV and calculate g|, H 
the moments as a function of P^ n ■ We- find a much smaller momentum dependence than the data 
indicate (see Figure Q). This suggests that the observed contribution from the high- mass states, 
i.e., the non-resonant X^f eso decays, is very important. 

In summary, for high lepton momenta, the moment measurements confirm earlier results by the 
CLEO Collaboration, and using A = 0.35 ±0.13 GeV we obtain Ai = -0.17 ±0.06 ±0.07 GeV 2 for 
P min = 1.5 GeV/ c. However, based on the OPE predictions for the P min dependence of the hadronic 
moments these values of the parameters fail to describe the moment measurements variation to lower 
lepton momenta. This inconsistency suggests that some of the underlying assumptions based on 
HQET and the OPE require further scrutiny. 
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